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Table 1: Ensemble details 
V aT x (GeV -1 ) crnifis) arrif(ud) 



am r , 



2.13 24 3 1.75(3) 0.04 0.03,0.02,0.01,0.005 0.0030(1) 
2.25 32 3 2.33(4) 0.03 0.008,0.006,0.004 0.00067(1) 



There are many reasons to study the light Baryons some of which are summarised by the late 
Nathan Isgur in "Why N*'s are important" [|lj], which I paraphrase here: 

• Nucleons are the stuff of which world is made. 

• They are the simplest system in which the quintessentially nonabelian character of QCD is 
manifest. There are N c quarks in proton because there are N c colours. 

• Baryons are sufficiently complex to reveal physics hidden from us in Mesons. 

Indeed, even the low lying spectrum is not completely understood. For example the nature of the 
Roper N(1440) J p = I resonance is a long standing puzzle. This state doesn't naturally fit into 
the quark models. In most models, the parity of the excitations of the nucleon are alternatively 
negative, then positive. The Roper has positive parity, and lies between Nucleon ground state 
N(939) J P = \ and the negative parity excitation N(1535) J p = \ . In the quark models, the 
positive parity excitation lies above the negative parity excitation and this has lead to speculation 
that the Roper is not an excitation of the nucleon but some other state. In principle lattice QCD can 
definitively resolve this issue simply by determining the spectrum of excited nucleons. However, 
excited states are difficult, requiring both large volumes and high statistics. 

The preliminary results presented in these proceedings are from the RBC and UKQCD 2+1 
flavor Domain Wall Fermion (DWF) data sets, with the Iwasaki Gauge action and the size of the 
fifth dimemsion, L$ = 16. The details are shown in table 111 the 24 3 data is described in |Q] and 
the 32 3 data in [||, ||]. DWF have lattice versions of the continuum QCD symmetries, in particular 
flavor, chiral and the lorentz (e.g. parity) symmetries. These symmetries come at an increased 
cost of simulation, but they protect matrix elements from mixing with other operators. This makes 
renormalisation of matrix elements simple 1 . So the target quantities for DWF calculations are 
matrix elements. Of course, one can use the same ensembles for spectrum calculations. Symmetry 
is also good for the spectrum, but what is really required is very high statistics, and what has been 
paid for with DWF is symmetry. In this sense, spectrum calculations are the poor relation in DWF 
calculations. Moreover, the smearings and sources tuned for the matrix element calculations are 
not optimal for the spectrum, and so each ensemble has different sources and numbers of sources. 

Shown in Figure [j] are the effective masses for the 32 3 data. As different sources were used 
to generate the data, the quality of the plateaux depends of the source used. In particular the 
mf(ud) = 0.006 was generated with a wall source and the data has a plateau which approaches 
from below. It is not clear the extent of the ground state saturation before noise swamps the signal. 
Certainly, the quality of the wall source data is inferior to the Gaussian source data. The results for 
the nucleon mass on the 24 3 data were presented at the previous lattice conference [pp. 

1 or at least, simpler. 
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Figure 1: Effection mass of the nucleon correlator for the 32 3 ensembles. The y-axis labels denote anif(ud). 
The upper and lower panels (circles) show correlators generated by the LHPC collaboration with Gaussian 
sources, the middle panel (squares) shows correlators generated with wall sources. The horizontal lines show 
the fit of an exponential to the correlator, and the span denotes the fit ranges. 

Ultimately the chiral extrapolation, to determine the spectrum at the physical quark masses, 
will need to be attempted, but its instructive to compare the nucleon masses of different simulations. 
This is done using the Edinburgh plot, and it is a good consistency test of a calculation. Shown in 
Figure ^ is the Edinburgh plot for all the 2+ 1 flavor DWF data. There are two things to note about 
this plot. Data from Lattice QCD now goes beyond the end of the quark model line. The data comes 
from different simulations and with changes to gauge action (DBW2 and Iwasaki), gauge coupling 
and hence lattice spacing, lattice volume and quark mass. At the level of statistical resolution the 
data lie on (apparently) a universal curve. One may conclude from this that the DWF calculations 
are internally consistent and these data are plotting the QCD curve. 

Many calculations now include sophisticated chiral extrapolations involving partially quenched 
data, with many different valence quark masses, including one on these data sets [||, Qj. Only uni- 
tary data, that is, the mass of the quarks in the sea is the same as the valence, is available for the 
baryon spectrum. Shown in Figure || are the nucleon masses. The simplest method for obtaining 
the continuum and chiral result is a linear extrapolation to the physical quark masses, and then 
"extrapolate" 2 the data to the continuum limit. However, more sophisticated chiral behavior can 
be examined. This is motivated by heavy baryon chiral perturbation theory |p], where the nucleon 
mass is given as a function of the pion mass squared. 
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"Naively drawing a straight line through two points. 
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Figure 2: The Edinburgh plot for all the DWF data generated on QCDOCs and at the Argonne Leadership 
Computing Facility, IBM BG/P (32 3 data). The solid line shows a quark model prediction 

This equation is used to suggest a form for a chiral fit for the nucleon mass at both lattice 
spacings, including quadratic lattice artifacts. 

3/2 ? 

M N = c + cim q + c 2 m q ' +c^a (2) 

This equation is then fitted to the five lightest data, am f = {0.004,0.006,0.008} from the 32 3 
ensembles and amf = {0.005,0.01} from the 24 3 ensembles simultaneously. This fit is also shown 
in Figure [| This a four parameter fit to five data and has an uncorrected # 2 /dof ~ 1.3, which is 
perhaps higher than one would like for an uncorrected fit. The nucleon mass then predicted is 

M N = 0.924(30)GeV (3) 

This error is statistical only as these are preliminary results, but does agree with the experimentally 
measured value. 

How sensible is this fit? Firstly, the curves plotted in Figure || are drawn using the fit param- 
eters with the lattice spacing set to its 24 3 value in red, 32 3 value in blue, and zero in black. The 
gradient of these curves is similar to the gradient of the straight line fits as they pass through the 
data. Whilst the gradient of the lattice spacing term has the opposite sign to the gradient of the 
continuum extrapololation shown in the inset, they are both relatively small, C3 ~ —0.17(5) for the 
chiral fit, ~ 0.57 for the linear fits. So lattice artefacts are small, again this can be seen in the plot 
as the data from different lattice spacings is close to each other. So, despite the large # 2 /dof, the 

3 /2 

fit reflects the data reasonable well. Moreover, the coefficent of the m q term, C2 can be related to 
gA once a fudge factor to convert from renormalised quark mass to pion mass squared is included. 
The value of gA obtained from the fit is then g^ ~ 2.2. This can be compared to the experimental 
value g^ xp ~ 1.27. This is clearly wrong, but critically, it is the right order of magnitude and the 
right sign. So, the fit seems reasonable given the data. 
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Figure 3: The nucleon mass versus renormalised quark mass in physical units. Blue (red) symbols denote 
32 3 (24 3 ) data. The straight lines show linear extrapolations at fixed lattice spacing. The curves show the 
chiral fit (equation (2)). The black symbol is the DWF prediction, and the green triangle is the experimental 
result. The inset shows the scaling behaviour of the linear extrapolation, and compares to experiment (green 
traingle) and the chiral fit (black diamond). 

Besides high statistics, large volumes are required for the baryon spectrum as these states are 
physically big, certainly compared to mesons. Finite size effects (FSE) are therefore an important 
source of systematic error. Several baryonic quantities have been measured and reported on these 
data sets. In ^jfy, a large FSE is reported for the 24 3 data for the axial charge of the nucleon. 
However, in [Q] no FSE was observed in the nucleon spectrum between the the 24 3 data and a 
smaller volume of 16 3 with all other quantities held fixed, for arrif = 0.01. The statistical precision 
of the data is around 1 — 2%. Subsequently, the lightest 24 3 ensembles have been extended, and it 
is these extended ensembles have been included in this work. The FSE is analyised for the heaviest 
quark mass in the chiral fit. However, it is instructive to compare the sizes of the lattice in terms 
of the Compton wavelength of the pion, m n L. For the 16 3 and 24 3 data used to estimate FSE, the 
m % L values are ~ 3.9 and ~ 5.6 respectively. For the lightest datum analysed in this work, the 32 3 
am/ = 0.004, the m K L value is ~ 4.1. This suggests that for the nucleon mass at least, FSE should 
be less than, say 1%. This is in accordance with "Lattice folklore" from other studies, that the FSE 
for the nucleon mass should be less than 1% for an m n L value of > 4. 

The nucleon operators used to construct the correlation functions are 

Q-i = (yCftVOV (4) 

n 2 = (ycy)v (5) 

On a lattice with anti-periodic boundary conditions, the forward propagating state of the correla- 
tion function constructed from the Hi operator projects onto the positive parity nucleon, and the 
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Figure 4: The effective masses of correlation functions constructed from the operators defined in equa- 
tions (4-5) with quark mass anif = 0.004. The label t + denotes the forward propagating state, and f_ denotes 
the time reversed backward moving state. 

backward propagtor which has the opposite parity projects onto the negative parity excitation, the 
N* state. The Q.2 operator has negative parity, and so the forward propagator of its correlation 
function projects onto a negative parity state and can be used as an additional estimate of the N*. 
The effective masses of these correlation functions are shown in Figure ||. The blue symbols show 
the effective masses for the negatice parity state, showing good agreement. Also shown (in red) 
is the backward propagator of Q.2 correlation function. What state does this project onto? It has 
positive parity, and is clearly not the nucleon. It is possible that this is the positive parity excitation 
of the nucleon. This state clearly lies above the negative parity excitation. 

Shown in Figure ^] is the light baryon spectrum for the 32 3 data. There are several limitations 
to this result, only a naive linear extraploation, only one lattice spacing, no estimate of FSE for the 
higher states, which could be potentially severe, and poor plateaux for the anif = 0.006 wall source 
data. However, it is interesting to compare to the physical spectrum, shown in lattice spacing units. 
The nucleon, delta and N* agree reasonably well with experiment. The higher excitations, the A* 
and the N** are more speculative, as they are more likely to be effected by the quality of plateaux 
and FSE. Considering the ordering of states only, if the backward moving Q.2 operator is indeed 
the positive parity excitation of the nucleon then this would exclude the N(1440) Roper resonance 
from being a nucleon. 

We have studied the light baryon spectrum for the 2+ 1 flavor DWF QCD ensembles and ex- 
plored some of the issues necessary to achieve the result. We present a preliminary result for 
the nucleon mass using a combined chiral/continuum fit which agrees with experiment. We also 
present a more speculative result for the excited spectrum, which taken at face value suggests that 
the Roper resonance is not a nucleon, but there several systematic uncertainties which are not suf- 
ficiently controlled to make this a concrete result. 
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Figure 5: Baryon masses in lattice units versus bare quark mass anif. 32 3 data only. Dotted vertical lines 
show the light quark mass and the chiral limit, labelled m m i and —am res respectively. Dashed horizontal lines 
show the physical spectrum, in lattice units. 
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